Investigation of Rock-Fluid Interactions in Carbonate Reservoirs with the Use of Digital Rock Physics

This thesis presents a pore-scale investigation into the mechanisms governing the interactions between water, oil and carbonate rocks using a series of core-flood experiments coupled with micro-CT imaging. Firstly, rock-fluid interactions leading to wettability alteration by seawater (injected in the secondary mode) and modified seawater (by addition or depletion of Ca2+ ion, injected in tertiary mode) are assessed using in-situ contact angle measurements. While the results have confirmed wettability alteration due to injection of the modified seawater (in tertiary mode), this change in wettability was found to be insufficient to fully explain the observed improvement of oil recovery. Secondly, a fluid-fluid interactions analysis methodology was developed and applied to the data generated as part of this thesis (Ca2+ ­enriched brine, Ca2+ depleted brine, SO42- enriched brine). The results show that emulsification is another mechanism responsible for enhancing oil recovery which is caused by reduction in interfacial tension. Data from multiple experiments (micro-CT imaging, stability analysis, surface tension measurements) also suggests that the amount of emulsion generated, and its stability depend on the type of ion added to seawater to modify the water chemistry. Specifically, Ca2+ depleted brine has shown to generate large amount of emulsion. As part of this thesis, a pore-scale direct flow simulation using the Volume-of-Fluid method was developed to model the mechanisms observed in the experimental studies, such as wettability alteration, viscous fingering and interfacial tension reduction. These mechanisms were grouped into two types, capillary-force affecting (contact angle and interfacial tension), and viscous-force affecting (viscous fingering reduction). When comparing the data between the two approaches (laboratory- and numerical experiments) the results show that capillary-force driven mechanisms were insufficient to predict the oil recovery due to development of flow stagnation zones in the sample. However, when the viscous forces were taken in account in estimation of sweep efficiency the displacement is significantly improved and the experimental data match well with that of numerical experiment. As a result, this thesis provides a first comprehensive pore-scale investigation of the physicochemical interactions between oil, water and carbonate rock related to the tertiary injections of modified seawater

Pore-GNN: A graph neural network-based framework for predicting flow properties of porous media from micro-CT images

This paper presents a hybrid deep learning framework that combines graph neural networks with convolutional neural networks to predict porous media properties. This approach capitalizes on the capabilities of pre-trained convolutional neural networks to extract n-dimensional feature vectors from processed three dimensional micro computed tomography porous media images obtained from seven different sandstone rock samples. Subsequently, two strategies for embedding the computed feature vectors into graphs were explored: extracting a single feature vector per sample (image) and treating each sample as a node in the training graph, and representing each sample as a graph by extracting a fixed number of feature vectors, which form the nodes of each training graph. Various types of graph convolutional layers were examined to evaluate the capabilities and limitations of spectral and spatial approaches. The dataset was divided into 70/20/10 for training, validation, and testing. The models were trained to predict the absolute permeability of porous media. Notably, the proposed architectures further reduce the selected objective loss function to values below 35 mD, with improvements in the coefficient of determination reaching 9%. Moreover, the generalizability of the networks was evaluated by testing their performance on unseen sandstone and carbonate rock samples that were not encountered during training. Finally, a sensitivity analysis is conducted to investigate the influence of various hyperparameters on the performance of the models. The findings highlight the potential of graph neural networks as promising deep learning-based alternatives for characterizing porous media properties. The proposed architectures efficiently predict the permeability, which is more than 500 times faster than that of numerical solvers.Document Type: Original articleCited as: Alzahrani, M. K., Shapoval, A., Chen, Z., Rahman, S. S. Pore-GNN: A graph neural network-based framework for predicting flow properties of porous media from micro-CT images. Advances in Geo-Energy Research, 2023, 10(1):39-55. https://doi.org/10.46690/ager.2023.10.0

Les droits disciplinaires des fonctions publiques : « unification », « harmonisation » ou « distanciation ». A propos de la loi du 26 avril 2016 relative à la déontologie et aux droits et obligations des fonctionnaires

The production of tt‾ , W+bb‾ and W+cc‾ is studied in the forward region of proton–proton collisions collected at a centre-of-mass energy of 8 TeV by the LHCb experiment, corresponding to an integrated luminosity of 1.98±0.02 fb−1 . The W bosons are reconstructed in the decays W→ℓν , where ℓ denotes muon or electron, while the b and c quarks are reconstructed as jets. All measured cross-sections are in agreement with next-to-leading-order Standard Model predictions.The production of $t\overline{t}$, $W+b\overline{b}$ and $W+c\overline{c}$ is studied in the forward region of proton-proton collisions collected at a centre-of-mass energy of 8 TeV by the LHCb experiment, corresponding to an integrated luminosity of 1.98 $\pm$ 0.02 \mbox{fb}^{-1}. The $W$ bosons are reconstructed in the decays $W\rightarrow\ell\nu$, where $\ell$ denotes muon or electron, while the $b$ and $c$ quarks are reconstructed as jets. All measured cross-sections are in agreement with next-to-leading-order Standard Model predictions

A study of CP violation in B-+/- -> DK +/- and B-+/- -> D pi(+/-) decays with D -> (KSK +/-)-K-0 pi(-/+) final states

A first study of CP violation in the decay modes $B^\pm\to [K^0_{\rm S} K^\pm \pi^\mp]_D h^\pm$ and $B^\pm\to [K^0_{\rm S} K^\mp \pi^\pm]_D h^\pm$, where $h$ labels a $K$ or $\pi$ meson and $D$ labels a $D^0$ or $\overline{D}^0$ meson, is performed. The analysis uses the LHCb data set collected in $pp$ collisions, corresponding to an integrated luminosity of 3 fb$^{-1}$. The analysis is sensitive to the CP-violating CKM phase $\gamma$ through seven observables: one charge asymmetry in each of the four modes and three ratios of the charge-integrated yields. The results are consistent with measurements of $\gamma$ using other decay modes

Study of the rare B-s(0) and B-0 decays into the pi(+) pi(-) mu(+) mu(-) final state

A search for the rare decays $B_s^0 \to \pi^+\pi^-\mu^+\mu^-$ and $B^0 \to \pi^+\pi^-\mu^+\mu^-$ is performed in a data set corresponding to an integrated luminosity of 3.0 fb$^{-1}$ collected by the LHCb detector in proton-proton collisions at centre-of-mass energies of 7 and 8 TeV. Decay candidates with pion pairs that have invariant mass in the range 0.5-1.3 GeV/$c^2$ and with muon pairs that do not originate from a resonance are considered. The first observation of the decay $B_s^0 \to \pi^+\pi^-\mu^+\mu^-$ and the first evidence of the decay $B^0 \to \pi^+\pi^-\mu^+\mu^-$ are obtained and the branching fractions are measured to be $\mathcal{B}(B_s^0 \to \pi^+\pi^-\mu^+\mu^-)=(8.6\pm 1.5\,({\rm stat}) \pm 0.7\,({\rm syst})\pm 0.7\,({\rm norm}))\times 10^{-8}$ and $\mathcal{B}(B^0 \to \pi^+\pi^-\mu^+\mu^-)=(2.11\pm 0.51\,({\rm stat}) \pm 0.15\,({\rm syst})\pm 0.16\,({\rm norm}) )\times 10^{-8}$, where the third uncertainty is due to the branching fraction of the decay $B^0\to J/\psi(\to \mu^+\mu^-)K^*(890)^0(\to K^+\pi^-)$, used as a normalisation.A search for the rare decays Bs0→π+π−μ+μ− and B0→π+π−μ+μ− is performed in a data set corresponding to an integrated luminosity of 3.0 fb−1 collected by the LHCb detector in proton–proton collisions at centre-of-mass energies of 7 and 8 TeV . Decay candidates with pion pairs that have invariant mass in the range 0.5–1.3 GeV/c2 and with muon pairs that do not originate from a resonance are considered. The first observation of the decay Bs0→π+π−μ+μ− and the first evidence of the decay B0→π+π−μ+μ− are obtained and the branching fractions, restricted to the dipion-mass range considered, are measured to be B(Bs0→π+π−μ+μ−)=(8.6±1.5 (stat)±0.7 (syst)±0.7(norm))×10−8 and B(B0→π+π−μ+μ−)=(2.11±0.51(stat)±0.15(syst)±0.16(norm))×10−8 , where the third uncertainty is due to the branching fraction of the decay B0→J/ψ(→μ+μ−)K⁎(892)0(→K+π−) , used as a normalisation.A search for the rare decays Bs0→π+π−μ+μ− and B0→π+π−μ+μ− is performed in a data set corresponding to an integrated luminosity of 3.0 fb−1 collected by the LHCb detector in proton–proton collisions at centre-of-mass energies of 7 and 8 TeV . Decay candidates with pion pairs that have invariant mass in the range 0.5–1.3 GeV/c2 and with muon pairs that do not originate from a resonance are considered. The first observation of the decay Bs0→π+π−μ+μ− and the first evidence of the decay B0→π+π−μ+μ− are obtained and the branching fractions, restricted to the dipion-mass range considered, are measured to be B(Bs0→π+π−μ+μ−)=(8.6±1.5 (stat)±0.7 (syst)±0.7(norm))×10−8 and B(B0→π+π−μ+μ−)=(2.11±0.51(stat)±0.15(syst)±0.16(norm))×10−8 , where the third uncertainty is due to the branching fraction of the decay B0→J/ψ(→μ+μ−)K⁎(892)0(→K+π−) , used as a normalisation.A search for the rare decays $B_s^0 \to \pi^+\pi^-\mu^+\mu^-$ and $B^0 \to \pi^+\pi^-\mu^+\mu^-$ is performed in a data set corresponding to an integrated luminosity of 3.0 fb$^{-1}$ collected by the LHCb detector in proton-proton collisions at centre-of-mass energies of 7 and 8 TeV. Decay candidates with pion pairs that have invariant mass in the range 0.5-1.3 GeV/$c^2$ and with muon pairs that do not originate from a resonance are considered. The first observation of the decay $B_s^0 \to \pi^+\pi^-\mu^+\mu^-$ and the first evidence of the decay $B^0 \to \pi^+\pi^-\mu^+\mu^-$ are obtained and the branching fractions, restricted to the dipion-mass range considered, are measured to be $\mathcal{B}(B_s^0 \to \pi^+\pi^-\mu^+\mu^-)=(8.6\pm 1.5\,({\rm stat}) \pm 0.7\,({\rm syst})\pm 0.7\,({\rm norm}))\times 10^{-8}$ and $\mathcal{B}(B^0 \to \pi^+\pi^-\mu^+\mu^-)=(2.11\pm 0.51\,({\rm stat}) \pm 0.15\,({\rm syst})\pm 0.16\,({\rm norm}) )\times 10^{-8}$, where the third uncertainty is due to the branching fraction of the decay $B^0\to J/\psi(\to \mu^+\mu^-)K^*(890)^0(\to K^+\pi^-)$, used as a normalisation

Search for the lepton flavour violating decay tau(-) -> mu(-)mu(+)mu(-)

A search for the lepton flavour violating decay $\tau^-\rightarrow\mu^-\mu^+\mu^-$ is performed with the LHCb experiment. The data sample corresponds to an integrated luminosity of 1.0 fb$^{−1}$ of proton-proton collisions at a centre-of-mass energy of 7 TeV and 2.0 fb$^{−1}$ at 8 TeV. No evidence is found for a signal, and a limit is set at 90% confidence level on the branching fraction, $\mathcal{B}(\tau^-\rightarrow\mu^-\mu^+\mu^-)<4.6\times10^{−8}$.A search for the lepton flavour violating decay τ$^{−}$ → μ$^{−}$ μ$^{+}$ μ$^{−}$ is performed with the LHCb experiment. The data sample corresponds to an integrated luminosity of 1.0 fb$^{−1}$ of proton-proton collisions at a centre-of-mass energy of 7 TeV and 2.0 fb$^{−1}$ at 8 TeV. No evidence is found for a signal, and a limit is set at 90% confidence level on the branching fraction, $\mathrm{\mathcal{B}}\left({\tau}^{-}\to {\mu}^{-}{\mu}^{+}{\mu}^{-}\right)<4.6\times {10}^{-8}$ .A search for the lepton flavour violating decay $\tau^-\to \mu^-\mu^+\mu^-$ is performed with the LHCb experiment. The data sample corresponds to an integrated luminosity of $1.0\mathrm{\,fb}^{-1}$ of proton-proton collisions at a centre-of-mass energy of $7\mathrm{\,Te\kern -0.1em V}$ and $2.0\mathrm{\,fb}^{-1}$ at $8\mathrm{\,Te\kern -0.1em V}$. No evidence is found for a signal, and a limit is set at $90\%$ confidence level on the branching fraction, $\mathcal{B}(\tau^-\to \mu^-\mu^+\mu^-) < 4.6 \times 10^{-8}$

Observation of the B0 → ρ0ρ0 decay from an amplitude analysis of B0 → (π+π−)(π+π−) decays

Proton–proton collision data recorded in 2011 and 2012 by the LHCb experiment, corresponding to an integrated luminosity of 3.0 fb−1 , are analysed to search for the charmless B0→ρ0ρ0 decay. More than 600 B0→(π+π−)(π+π−) signal decays are selected and used to perform an amplitude analysis, under the assumption of no CP violation in the decay, from which the B0→ρ0ρ0 decay is observed for the first time with 7.1 standard deviations significance. The fraction of B0→ρ0ρ0 decays yielding a longitudinally polarised final state is measured to be fL=0.745−0.058+0.048(stat)±0.034(syst) . The B0→ρ0ρ0 branching fraction, using the B0→ϕK⁎(892)0 decay as reference, is also reported as B(B0→ρ0ρ0)=(0.94±0.17(stat)±0.09(syst)±0.06(BF))×10−6

Measurement of the CP-violating phase β\beta in B0→J/ψπ+π−B^0\rightarrow J/\psi \pi^+\pi^- decays and limits on penguin effects

Time-dependent CP violation is measured in the $B^0\rightarrow J/\psi\pi^+\pi^-$ channel for each $\pi^+\pi^-$ resonant final state using data collected with an integrated luminosity of 3.0 fb$^{-1}$ in $pp$ collisions using the LHCb detector. The final state with the largest rate, $J/\psi\rho^0(770)$, is used to measure the CP-violating angle $2\beta^{\rm eff}$ to be $(41.7\pm 9.6_{-6.3}^{+2.8})^{\circ}$. This result can be used to limit the size of penguin amplitude contributions to CP violation measurements in, for example, $B_s^0\rightarrow J/\psi\phi$ decays. Assuming approximate SU(3) flavour symmetry and neglecting higher order diagrams, the shift in the CP-violating phase $\phi_s$ is limited to be within the interval [$-1.05^\circ$, +$1.18^\circ$] at 95% confidence level. Changes to the limit due to SU(3) symmetry breaking effects are also discussed.Time-dependent CP violation is measured in the B(−−−)0→J/ψπ+π− channel for each π+π− resonant final state using data collected with an integrated luminosity of 3.0 fb −1 in pp collisions using the LHCb detector. The final state with the largest rate, J/ψρ0(770) , is used to measure the CP -violating angle 2βeff to be (41.7±9.6−6.3+2.8)° . This result can be used to limit the size of penguin amplitude contributions to CP violation measurements in, for example, B(−−−)s0→J/ψϕ decays. Assuming approximate SU(3) flavour symmetry and neglecting higher order diagrams, the shift in the CP -violating phase ϕs is limited to be within the interval [ −1.05°,+1.18° ] at 95% confidence level. Changes to the limit due to SU(3) symmetry breaking effects are also discussed.Time-dependent CP violation is measured in the $B^0\rightarrow J/\psi\pi^+\pi^-$ channel for each $\pi^+\pi^-$ resonant final state using data collected with an integrated luminosity of 3.0 fb$^{-1}$ in $pp$ collisions using the LHCb detector. The final state with the largest rate, $J/\psi\rho^0(770)$, is used to measure the CP-violating angle $2\beta^{\rm eff}$ to be $(41.7\pm 9.6_{-6.3}^{+2.8})^{\circ}$. This result can be used to limit the size of penguin amplitude contributions to CP violation measurements in, for example, $B_s^0\rightarrow J/\psi\phi$ decays. Assuming approximate SU(3) flavour symmetry and neglecting higher order diagrams, the shift in the CP-violating phase $\phi_s$ is limited to be within the interval [$-1.05^\circ$, +$1.18^\circ$] at 95% confidence level. Changes to the limit due to SU(3) symmetry breaking effects are also discussed.Time-dependent CP violation is measured in the B0→J/ψπ+π− channel for each π+π− resonant final state using data collected with an integrated luminosity of 3.0 fb −1 in pp collisions using the LHCb detector. The final state with the largest rate, J/ψρ0(770) , is used to measure the CP -violating angle 2βeff to be (41.7±9.6−6.3+2.8)° . This result can be used to limit the size of penguin amplitude contributions to CP violation measurements in, for example, Bs0→J/ψϕ decays. Assuming approximate SU(3) flavour symmetry and neglecting higher order diagrams, the shift in the CP -violating phase ϕs is limited to be within the interval [ −1.05°,+1.18° ] at 95% confidence level. Changes to the limit due to SU(3) symmetry breaking effects are also discussed